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Ph. Boucaud'^, F. De Soto*^, A. Donini^ J.P. Lcroy^, A. Le Yaouanc'^, J. Micheli'', H. Moutarde'^, O. 
Pene^, J. Rodn'guez-Quintero" 

'^Laboratoire de Physique Thcorique, Universite de Paris XI, Batimcnt 210, 91405 Orsay-Cedex, France. 

'^Dpto. de Ffsica Atomica Molecular y Nuclear, Universidad de Sevilla. Avda. Rcina Mercedes s/n, 
41012 Sevilla, Spain. 

'^I.N.F.N., Roma I and Dip. Fisica, Universita di Roma "La Sapienza", P. le A. Moro 2, 00185 Rome, 
Italy 

'^Centre de Physique Theorique Ecole Polytechnique, 91128 Palaiseau Cedex, France 
°Dpto. de F.A., E.P.S. La Rabida, Universidad de Huelva, 21819 Palos de la Fra., Spain. 

Gluon two- and three-point Green Functions computed in Landau gauge from the lattice show the existence 
of power corrections to the purely perturbative expressions, that can be explained through an Operator Product 
Expansion as the influence of a non gauge invariant mass dimension two condensate. The relationship of this 
condensate with topological properties of QCD, namely instantons, will be studied, giving a first estimate of the 
contribution of instantons to this condensate based in the direct lattice measure, after a cooling process, of the 
instanton liquid properties. 



1. QCD coupling constant, O.P.E. and 
(^■^) condensate. 

Lattice calculations of the QCD coupling con- 
stant and gluon propagator based in the Green 
Functions Method [Q, suggest the necessity to 
add power corrections to the purely perturbative 
expressions to correctly describe their running . 
An Operator Product Expansion (O.P.E.) analy- 
sis of the Green functions in Landau gauge|^ re- 
lates this power corrections to the existence of 
a non-perturbative (^^) condensate Q, through 
expressions: 

O.P.E. iP ) ^ '^Pert.iP ) +C ^2 ' 

^In the lattice we will work in the minimum Landau 
gauge, di^Afi = 0, so all gauge dependent quantities, will 
be expressed in this particular gauge. 



where perturbative expressions are developed at 
three loops, and the functions c and c' include 
the Wilson coefficient of the expansion, and the 
anomalous dimension of the condensate at leading 
logarithm. 

By performing a combined fit of lattice re- 
sults to expressions in (|l|), in two different MOM 
schemes, a value of Aj^ is extracted, in fairly 
good agreement with the one obtained by the AL- 
PHA collaboration by a completely different 
method. A value of the condensate comes out 
from the analysis. The physical meaning of this 
condensate is still an open question, and a lot 
of work is being devoted to its study during last 
years, for example, in relation to confinement . 
The aim in this work will be to study the possi- 
ble semiclassical contribution to this condensate 
coming from instantons, and whether they might 
explain the presence of power corrections in Green 
Functions. 
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2. The role of instantons. 

Instantons have been extensively studied as a 
possible description of the QCD vacuum (See Q 
for a general overview), and so as a major source 
of QCD properties at low energies. In relation 
with the aim of this work, an ensemble of non- 
interacting instantons (I) and antiinstantons (T) 
in Landau gauge would give a contribution to the 
(A^) condensate; 
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- I d''xA;ix)A;{x) = Un'p'n, (2) 



where A'^{x) is the standard 't Hooft Polyakov 
instanton gauge field [^, p the average radius, 
and n = ^'^^i the density. 

If we accept the phenomenological values as- 
signed to n and p by the Instanton Liquid Model 

(n ~ 0.5fm~^ and p ~ 1/3/to), the instan- 
tonic contribution will be (A'^), , ^ 0.5GeV^. 

\ I Inst. 

We will perform, however, our own analysis, thus 
testing the latter approach. 

2.1. Cooling. 

In principle, a direct measure of A^ in the lat- 
tice should be possible, but the presence of the 
UV divergent part is hardly separable from the 
soft, instantonic one. The other possibility is to 
perform a cooling procedure, that will allow us 
to compute the number and size of instantons, 
giving an indirect measure of the A'^ through (0). 

We will use the traditional cooling method |||], 
even if it introduces a number of known biases, as 

1 — I annihilation, and a modification of instanton 
sizes and lattice spacing. The approach proposed 
here is to compute instanton properties for dif- 
ferent number of cooling sweeps, and extrapolate 
back to the thermalised situation in order to re- 
cover their physical meaning^. 

2.2. Shape Recognition. 

Instantons will be localised in cooled lattices 
via a geometrical method (Described in pO[.) 



^ The use of improved cooling methods, as the one devel- 
oped in [b|, could improve this approach, as radii evolution 
is minimised, but I — I annihilation is unavoidable, so the 
extrapolation will be anyway necessary. 



that accepts a topological charge lump as an in- 
stanton when the ratio of the integral over a given 
fraction of the topological charge at the maxi- 
mum, a, and its theoretical counterpart, e, is 1, 
for a range of values of a: 



^^,\Qpm>^d'^xQp{x 

I IQp(0)l - 

1 - 3ai/2 + 2a3/4 



(3) 



Once the lump has been identified as an instan- 
ton, the radius will be computed from the size of 
the cluster where the integral has been developed. 

2.3. A naive model of annihilation. 

With the method outlined above, we compute 
the density and size of instantons in a lattice, for 
different numbers of cooling sweeps, tIc, obtain- 
ing values with a strong dependence on rtc (See 
figure), that avoids to obtain any physical infor- 
mation at fixed Uc- 

As a first approach to the understanding of this 
evolution, we will make a simple model, where 
instantons annihilate with antiinstantons (Being 
so AA = Nj — Nj a constant) proportionally to 
their packing ratio, and to the number of antiin- 
stantons, so that the equation for the evolution 
of A = A/ + Aj is: 



If we assume p{nc) — cte, the solution of Eq. 



(4) 



will give N{nc) ^ j^^^^, the expression used in 
pot , as a first order approach, but our coohng 
procedure modifies instanton's size (See figure), 
in a way thah we phenomenologically parametrise 



p{nc) = p(0)(l + aln(l + ?ic)) 



C5) 



We will include (||) in equation (^), with p{Q) the 
extrapolated radius at the thermalised situation 
and a a constant to determine. 

After performing a combined fit of our lattice 
results to the expressions (^) and the one coming 
from the integration of (|^), we can fix the initial 
values of the density, n(0) and the radius p{0), 
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and the two constants that govern the evolution, 
A and a. 
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Figure 1. Results of the combined fit for the instan- 
ton density and radius as a function of the number of 
cooling sweeps for a 24* lattice at (3=6.0. 

3. Results and conclussion. 

The result of the combined fit gives a value 
of the instantonic contribution to (^/„s) ~ 
0.4GeV^, however the result of the extrapolation 
is highly dependent on the value of p, which due 
to the logarithmic behaviour is hardly reliable. 
We therefore prefer the value at the maximum, 
1.12(11) GeV^, as a crude estimation of (^j„g). 

This semiclassical evaluation of (^^), which 
does not run with the scale, is difficult to re- 
late to that appearing in the O.P.E. expansion, 
which does depend on the renormalisation scheme 



and scale. The typical scale of instantons is 
p^^ ~ 0.7GeV. Unluckly it is not possible to run 
the (j4q ) to such a low energy, where per- 
tubative QCD is not valid. The lowest reacheble 
energy scale is 2.6GeV [|,|lO|; 

(^2, p^.(2.6GeV)) = 1.4(3)(3)GeV2, (6) 

the first error coming from the OPE determina- 
tion of the condensate renormalised at 10 GeV, 
and the second from higher orders in the running. 

Keeping in mind the level of uncertainty of 
these calculations, we can nevertheless claim a 
rather encouraging agreement between the in- 
stantonic contribution to the condensate and the 
one computed from the running of the Green 
Functions. 
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